Nonmonotonic gap in the coexisting antiferromagnetic and superconducting state for 

electron-doped cuprate superconductors 
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We argue that the experimentaUy observed nonmonotonic gap in electron-doped cuprates at 
optimal doping is the lowest quasiparticle excitation energy in the coexisting antiferromagnetic (AF) 
and superconducting (SC) state. The idea is implemented by studying the coexistence of AF and 
SC orders with the t-t'-t"-J model. Although the pairing gap itself is assumed to be the simplest 
d wave which is monotonic, we have found that the quasiparticle excitation gap in the coexisting 
state is nonmonotonic, with the maxima around the hot spots where the Fermi surface is missing 
due to the AF gap. Within the same framework of the coexisting state the spectral function is also 
calculated at optimal doping. The obtained results are all consistent with experiments. 
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The pairing symmetry is a central issue for under- 
standing the notable superconducting (SC) properties 
in cuprate superconductors. While for hole-doped com- 
pounds the pairing symmetry is generally accepted to 
be d wave, it is still under debate for electron-doped (e- 
doped) ones. Originally it was thought to be s wave,^ 
but later suggested to be d wave as in hole-doped ma- 
terials, by a number of experimental measurements such 
as phase sensitive,^ angle-resolved photoemission spec- 
troscopy (ARPES)r' penetration depthj"* etc. On the 
other hand, recent experimental evidence disfavoring the 



typical d wav( 



■ 5.6.7 



or supporting a transition from d 



to s wave with increasing doping was also found. In 
this regard, it is of great interest to notice the measure- 
ments by Raman scattering^ and ARPESS on e-doped 
compounds Nd2_xCe2:Cu04 (NCCO) at x = 0.15 and 
Pri_2:LaCe2:Cu04 (PLCCO) at x = 0.11, respectively. 
They are both on optimally doped samples and have 
reached the same conclusion, i.e., the observed gap does 
not fit the simplest commonly assumed d-wave function 
coskx — cosfcy, but exhibits a nonmonotonic behavior 
with the gap maxima locating midway between the Bril- 
louin zone (BZ) boundaries and the zone diagonals. 

So far, the theoretical explanations for the so called 
nonmonotonic d-wave gapSiL^iiSiii are limited to extend- 
ing the SC gap out of the simplest d wave, e.g., by inclu- 
sion of high-order d-wave harmonics. 

Here an alternative idea is proposed to explain the 
observed nonmonotonic gap. The new idea highlights 
the antiferromagnetic (AF) order in the SC phase, which 
is motivated by the following facts. First, in e-doped 
cuprates the AF order is robust to survive a broad dop- 
ing range, which coexists with the SC order around the 
AF/SC phase boundary and even at optimal doping. ^^'^^ 
Actually, for the SC samples of NCCO at a; = 0.15, the 
Neel temperature Tn is usually much higher than the SC 
transition temperature Tci^^i^^ Second, recent ARPES 
measurements have revealed the intriguing doping evolu- 
tion of the Fermi surface (FS) in NCCO>i& It was found 
that at optimal doping the FS consists of two inequiva- 



lent pockets around (tt, 0) and {tt/2,tt/2). This suggests 
the two-band modeling to study the SC properties in e- 
doped cuprates, as recently phenomenologically adopted 
by Luo and Xiang to calculate the superfluid densitypi^ 
Theoretically the presence of two inequivalent FS pock- 
ets (and two bands) has been attributed to the band 
folding due to the AF order pi2*i2iSfliSi From all of these, 
we argue that the experimentally observed gap at opti- 
mal doping is the lowest quasiparticle excitation energy 
in the coexisting AF and SC state. Due to the AF gap 
the FS is segmented and partly missing as seen from the 
ARPES, thus the lowest quasiparticle excitation energy 
is not purely the SC gap. Then, even if the SC gap itself 
is assumed to be the typical d wave which is monotonic, 
the quasiparticle excitation gap may be nonmonotonic. 

With the above idea in mind, in this paper we inves- 
tigate the coexistence of the AF and SC orders in e- 
doped cuprates, based on the t-t'-t"-J model. Although 
it is currently not reconciled whether the strong-coupling 
model is suitable for e-doped cuprates (Ref. 20] and ref- 
erences therein), it provides a concise realization of the 
idea which is itself independent of the models. As usual, 
the SC pairing from the J term is assumed to have the 
typical d-wave symmetry. Despite this, we have found 
that the quasi-particle excitations in the coexisting state 
have a nonmonotonic gap, with the maxima around the 
hot spots where the FS is missing. This is consistent with 
the experimental observationn^ Also, the spectral func- 
tion at optimal doping for NCCO is calculated for the 
first time within the coexisting AF and SC state. The 
obtained result is in agreement with the ARPES data>i& 

We begin with the t-t'-t"-J model Hamiltonian on a 
square lattice 



(ij)-2(T 



-t" E + H.c.) -f jY,{S, ■ S, ~ -n,n,) 
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where (), ()2, Os represent the nearest neighbor (n.n.), 
second n.n., and third n.n. sites, respectively, and the 
rest of the notation is standard. No double occupancy is 
allowed in the model, which implies that the Hamiltonian 
is already the transformed version after particle-hole 
transformation. Thus electron-doping already becomes 
hole-doping, and correspondingly one has t < 0, t' > 
and t" < 0. Throughout the work \t\ is taken as the 
energy unit. 

The slave-boson transformation is used to treat the 
Hamiltonian, i.e., c^cr — blfi„ with 6^: bosonic holon 
operator, fi^: fcrmionic spinon operator and the con- 
straint blbi + fiaficr = 1 at cach site. Then boson 
condensation is assumed, i.e., (bi) = {b\) = ^/x {x: dop- 
ing concentration) at low temperatures we are interested 
in. Subsequently the left spinon Hamiltonian is decou- 
pled with definition of the following mean-field param- 
eters: the uniform bond order (fj^-fja) = Xa = Xi the 
AF order {fj-^fi^ - fljii)/'^^ = (-l)*m, and the spinon 
pairing order [or resonating-valence-bond (RVB) order] 
ifi]'fjl ~ filfj"\) ~ ^iji where the standard d-wave sym- 
metry is assumed, i.e., A^j = A (—A) for bond {ij) along 
X (y) direction. Note that the RVB order is directly as- 
sociated with the SC order at finite doping as discussed 
below. We mention that the technical treatment is simi- 
lar to that used early by Inui et ali^ Their work, however, 
does not include t' and t" which are essential to exhibit 
the electron-hole asymmetry and crucial to all the re- 
sults discovered here. In consideration of two sublattices 
D and E (with corresponding operators d and e, respec- 
tively), the Hamiltonian after the above treatment, 
can be expanded as follows in momentum space (up to 
irrelevant constants) 



And the free energy (per site) is written down 

F/N = -{2T/N) J2 ln[2cosh(ii;,72T)] 

—fix 



(5) 



For each given doping x and temperature T, the order 
parameters x, to, A and the chemical potential /i are 
determined by minimizing the free energy. 
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FIG. 1: Energy bands at doping x = 0.15: (a) quasiparticle 
bands in the presence of both AF and SC orders (with 



A = 0.0916, m = 0.144, 



X 



(b) hypothetical AF bands 



(2) 



where 
Afc 



= {— 2tx — J x) {cos kx + cos ky), 
= —At'x cos kx cos ky — 2t"a;(cos 2kx 
~ A(cosfca; 



^y 

COS ky), 



COS 2k„ 



fi is the renormalized chemical potential and N is the 
total number of lattice sites. Note that the wavevector k 
is restricted to the magnetic Brillouin zone (MBZ). 

The Hamiltonian 10) can be diagonalized, giving rise 
to the quasiparticle energy bands i-E^, with 



4± 



4J2 



TO^ 



(3) 
(4) 



0.189, M = -0.00385) and 
fj, in the absence of pairing 
(with the same m, x ^-nd /i as above). In each panel the 
solid (dashed) line represents the band '4-(— )'. The model 
parameters are: t' = 0.32, t" = —t'/2, J = 0.3 (in units of 
\t\)- 

In the following we are limited to T ~ 0. And we are 
mainly interested in optimal doping, i.e., x — 0.15 for 
NCCO. It is expected at this doping to obtain a coexist- 
ing AF and SC state which has nodal quasiparticle exci- 
tations. This can be easily realized with choice of typical 
model parameters: t' = 0.32, t" ~ —0.16 and J — 0.3, 
leading to A = 0.0916, m = 0.144, x ^ -0.189, /i = 
—0.00385 at a; = 0.15. The quasiparticle energy bands 
, defined in the MBZ, are extensively plotted in the 
original BZ in Fig.^a). The zero energy excitations, con- 
tributed by the band, are present around (7r/2,7r/2) 
(see the inset). This will lead to the 'V-shaped' density of 
states at very low energy. In order to clearly understand 
the formation of the quasiparticle bands, we have artifi- 
cially plotted the energy bands — fJ- with the same pa- 
rameters TO, X Slid fj, as above, which can be regarded as 
the AF bands in the absence of pairing. Both bands are 



3 



crossed by the Fermi level, around (7r,0) and (7r/2,7r/2), 
respectively. The corresponding FS is plotted in Fig. |21 
in the first quadrant by the red solid and dashed lines. 
It consists of three pieces, which are separated due to 
the AF gap, in contrast to the continuous FS around 
(tt, tt) for free particles as shown by the thick black line 
in Fig. 121 When the pairing is considered, the RVB gap 
will be further opened, leading to the quasiparticle bands 
in Fig.^a). 
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FIG. 2: Fermi surfaces in the first quadrant. The red lines 
give the FS for the AF bands in Fig.^b), of which the solid 
line is contributed by band '+' and the dashed ones by ' — '. 
[The FS for the ' — ' band has been moved outside of the MBZ 
by wavevector (7r,7r).] For comparison, the thick black line 
around (7r,7r) shows the FS for free particles (without AF 
order). Its crosses with the MBZ boundary are so called hot 
spots, as indicated by the black dots. 



(/i in Fig. 121 the momentum components and ky satisfy 
the relation tan0 — (jr — k^) / (tt — ky) . In each direction, 
the energy minima Ag^i and Ag^2 for bands E'^ and E^^ , 
respectively, can be found out, which are shown in the 
upper inset of Fig. |3| The final excitation gap in the 
direction cj) should be the smaller one of Ag^i and Ag^2, 
as given by the main panel of Fig. O It is clear that 
Ag shows a nonmonotonic behavior, with the gap max- 
imum at (/)' = 18°, around one of the hot spots shown 
in Fig. 121 This is qualitatively the same as the experi- 
mental results)^ see the lower inset of Fig. 13 where the 
data from Ref. are re-plotted. The maximum position 
(j)' = 18° is close to 4>' — 15° or 20° observed in optimally 
doped NCCO^ or PLCCO^ respectively. 

The nonmonotonic gap can be explained as follows 
based on the FS (red fines) in Fig.|21 At = 45° (0' = 0), 
i.e., in tfie zone diagonal direction, the FS contributed by 
the band is crossed and at the same time the d-wave 
pairing gap vanishes, thus the excitation energy is zero. 
When (j) slightly deviates from the nodal direction, the 
FS is still crossed and thus the excitation gap is just the 
pure pairing gap for the band. Once (f) decreases so 
that the direction points to the region where the FS is 
missing, the excitation gap is combined from the nonzero 
term — fJ. and the pairing gap [see Eq. ©]. Due to the 
extra energy — fi\ caused by the AF gap, Ag increases 
rapidly and reaches the maximum when A^.i = Ag^2- 
Then Ag transfers from Ag^i to Ag^2- Still, the current 
Ag^2 is combined from the term — fi and the pairing 
gap. When continues to decrease, the FS contributed 
by the ' — ' band is crossed. And Ag becomes again the 
pure pairing gap for the '— ' band, which is also d wave, 
and persists tifi (/) = ((/>' = 45°). 
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FIG. 4: Spectral intensity at doping a; = 0.15 in the coexist- 
ing AF and SC state, obtained by integrating A{k, ui) times 
Fermi function over an energy interval [—40, 20] meV. The 
Lorentzian broadening is 0.03|t|. 



FIG. 3: Quasiparticle excitation gap Ag (in unit of |t|) as 
a function of 0'. Upper inset: energy minima of bands E'^ 
(solid line) and (dashed line) in each direction 0'. Lower 
inset: a sketch of the experimental gap for NCCO at a:: = 0.15, 
extracted from Ref. 0. 

From the energy bands E^ , the excitation gap Ag can 
be extracted. For a given direction as indicated by angle 



Within the same framework of the coexisting state, we 
further calculate the single particle spectral function, in 
order to compare with the ARPES result on NCCO at 
optimal doping x — O.lbJ^ The similar calculation at 
optimal doping was previously undertaken, but limited 
to the pure AF stateiiSiiSiSS To match the experimental 
situation at a; = 0.15 and low temperature T < Tj.^ the 
coexisting state including the SC order is necessary. For 
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primitive physical electrons, which are represented by op- 
erators Ckcr [equal to cj^^ in Hamiltonian ^] , the Matsub- 
ara Green's function is G„{k,T) = —{TrCka{T)cl^{0)) — 
~x{TT-fl^{T)fka{0))- The Fourier transformation of G is 
eventually derived, irrelevant of spin ct, as follows 



G{k,iujn) 



(l-Pk) 



l-r: 



El 



iujn - E, 



+ (I+Pk) 



1 



(6) 



with the fcrmionic Matsubara frequencies ujn and the co- 
efhcients 



Pk 



el + 4j2r 



ri = i^i-f^)/E^ 



The spectral function is A{k,uj) = 
— (l/7r)ImG(fc, iw„)|j;^^_,j^^jo+. For direct compari- 
son with experiment, we have calculated the integra- 
tion JdLjA{k,Uj)/{l + over an energy interval 
[-0.12,0.06]|t| ~ [-40,20] meV and drawn its density 
plot in Fig. ^ Strong intensity is concentrated on several 
separate regions, i.e., the patch around (7r/2, 7r/2), a half 
of the pocket around (tt, 0) and equivalently that around 
(0, tt). All of them, when looked together, are close to a 
large curve around (tt, tt). This is well consistent with 
the ARPES result (Fig. 3 of Ref. O for x = 0.15). 

Compared to the corresponding result (not shown) cal- 
culated in the pure AF state embodied in Fig.^b), the 
spectral intensity obtained here in the coexisting state is 
qualitatively similar, but quantitatively weaker mainly 
in the region around (tt, 0) [and (0,7r)]. This is un- 
derstandable in view that the d-wave SC gap is actu- 
ally small around (7r/2,7r/2) and only becomes notable 
around (7r,0). Thus the current result has justified the 
part reasonability of the previous calculations at optimal 



doping which were based on the pure AF stateii2ii2i2£ 

On the other hand, it can be imagined that if the pair- 
ing were much stronger (i.e., the order parameter A were 
much bigger) the spectral intensity would vanish around 
(tt, 0), but remain around (7r/2,7r/2). 

A few remarks are in order. First, the simplest d-wave 
pairing is assumed here. It might be true that the pairing 
itself has a more complex symmetry, which is, however, 
not the emphasis of our idea. Second, the mean-field 
treatment tends to overestimate the AF order and per- 
haps simultaneously the SC order. Thus smaller order 
parameters m and A are expected under consideration 
of the quantum fluctuations, leading to a reduced ex- 
citation gap as shown in Fig. 13 Finally, the spatially 
uniform coexistence of both AF and SC orders is implied 
in our calculation. Experimentally it is not conclusive^— 
whether the two orders are uniformly coexistent or not. 
In the latter case, further study is needed to extract the 
excitations in the inhomogeneous phase. 

In conclusion, we have studied the coexistence of the 
AF and SC orders within the t-t'-t"-J model for e-doped 
cuprates. Although the pairing gap is assumed to be the 
standard d wave which is monotonic, we have found that 
the quasiparticle excitation gap in the coexisting state 
is nonmonotonic, with the maxima around the hot spots 
where the FS is missing due to the AF gap. This provides 
an explanation for the nonmonotonic gap observed at 
optimal doping. The calculated spectral intensity in the 
coexisting state is also consistent with ARPES result. 
Our work suggests that the e-doped cuprates at optimal 
doping could be understood within a unified theory based 
on the AF and SC coexistence. 
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